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I. INTRODUCTION 


Hydrodynamic equations are conventionally expressed 


Abstract— This study serves as an extension of prior research focusing 
on weighting coefficients within the context of weighted Taylor series. 
The primary objective is to determine the weighting coefficients’ values 
in the weighted Taylor series for the purpose of modeling water waves 
based on velocity potential. Utilizing the weighted Taylor series, we 
derive both the weighted continuity equation and the weighted Laplace 
equation. The latter is addressed using the variable separation method 
over a sloping bottom, leading to the formulation of the velocity 
potential equation, wave constant equations, and energy conservation 
equations. Within the wave constant equations, a breaking equation is 
incorporated. Leveraging both the breaking equation and the energy 
conservation equations, breaking indexes equations are formulated. 
These equations encompass breaker length, breaker depth, and breaker 
height indexes, with weighting coefficients prominently featured. 
Calibrating the results of the breaking indexes equations against 
findings from earlier studies provides suitable values for the weighting 
coefficients. Additionally, this research introduces a shoaling-breaking 
model and a refraction-diffraction model to explore the phenomena of 
shoaling-breaking within the solution of the weighted Laplace equation. 


order differential terms, ensuring their disappearance 


without merely being eliminated or overlooked. 


through Taylor series, typically truncated to first order, 
under the assumption that higher-order terms, such as 
second order and beyond, become negligible at very small 
intervals. While not inherently incorrect, this truncation 
approach sacrifices the representation of the function's 
characteristics encapsulated in the higher-order terms, 
leading to an imprecise approximation. 


To address this limitation, Hutahaean (2021,2022,2023a) 
proposes the utilization of weighted Taylor series. This 
involves truncating the Taylor series to the first order, while 
compensating for the of higher-order 
contributions through the incorporation of weighting 
coefficients. work (2023a), 
systematic application of Taylor series truncation is 
introduced, employing the Central Difference Method. This 
method allows for the systematic removal of even higher- 


omission 


In Hutahaean's a more 
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Furthermore, the study explores the extraction of 3rd order 
odd differential term contributions, with the additional 
inclusion of 5th order terms. 


The use of the weighted Taylor series is characterized by 
the absence of truncation errors. Consequently, when 
applied to the continuity equation, the weighted Taylor 
series formulation minimizes or eliminates truncation 
errors. Similarly, the weighted Laplace equation, derived 
from the weighted continuity equation, maintains this 
heightened accuracy. Moreover, equations pertaining to 
various wave mechanics, formulated using the weighted 
Laplace equations, incorporate weighting coefficients, 
thereby enhancing their accuracy and reliability. 
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Il. WEIGHTED TAYLOR SERIES 


Taylor series (Arden, Bruce W. and Astill Kenneth N. 
,1970) for a function f = f (x,t) where x is the horizontal 
axis and t is time 


(x + 6x,t + 6t) = f( grat Lys af y 
PETA = f(x, at ax 
ôt? d? da? ôx? d? 
Ss stax f E 
2! dt? dtdx 2! dx? 
SEPRE (1) 


This equation can be expressed as, 


f(x + 6x,t + 6t) = f(x,t) + 


std a af 
(V+ patoret ott 


2! dt dx dt 
ôx d df 

(ie a 

( i 2! ax ) d 


In this paper, x represents the horizontal axis, z represents 
the vertical axis, while t represents time. 


The complete Taylor series, there are contributions from 
high-order differential terms in the first-order differential 
term. The contributions of these high-order differential 
terms can be represented by a coefficient, allowing for the 
elimination of terms with high differentials. 


df df 
f(x + 6x,t + dt) = f(x,t) + Ye2ôt mr + YxÔX a 


Yt,2 and y, are referred to as weighting coefficients, while 
(2) is called the weighted Taylor series. The weighted 
Taylor series for the function f (x, z, t) is as follows: 


d 
f(x + 6x,z + 6z,t + ôt) = f(x,z,t) treast 


af af 
+yxôx aT Vz0Z ae (3) 


There is no difference in the values of y, in the function 
f(x,t) compared to yx in the function f(x, z,t). The 
baseline values for these weighting coefficients are yp,2 = 
2, Yt3 = 3, Yx = 1 and y; = 1. More precise values of the 
weighting coefficients are presented in Table (1). The 
calculation method can be found in Hutahaean (2023a), 
where the computation of the weighting coefficients in 
Table (1) is achieved with greater accuracy compared to 
Hutahaean (2023a). In these weighting coefficients, all 
even-order differential terms in the Taylor series are 
represented, while the absorbed odd-order high differential 
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terms are of order 3 and 5. Higher-order differentials are 
truncated with the reduction of intervals, dt, 6x and 6z. The 
values of the weighting coefficients are presented in Table 


(1). 
Table 1 Weighting coefficient values 


E Yt,2 Yt,3 Yx Yz 

0.03 | 1.99812 | 3.05087 | 0.98899 | 1.11999 
0.032 | 1.99786 | 3.05877 | 0.98746 1.139 
0.034 | 1.99758 | 3.06738 | 0.98583 | 1.15986 


0.036 | 1.99728 | 3.07674 0.9841 | 1.1827 
0.038 | 1.99695 | 3.08689 | 0.98227 | 1.20765 

0.04 | 1.99662 | 3.09786 | 0.98034 | 1.23485 
0.042 | 1.99626 3.1097 | 0.97831 | 1.26448 
0.044 | 1.99588 | 3.12244 | 0.97617 | 1.29669 


0.046 | 1.99548 | 3.13614 | 0.97393 | 1.3317 


€ is referred to as the optimization coefficient, where a 
larger € corresponds to a larger interval size and 
consequently, larger values for the absorbed third and fifth- 
order differential terms in the weighting coefficients. This 
research aims to obtain appropriate values for the weighting 
coefficients expressed in terms of £. The appropriateness is 
assessed concerning the breaking parameter generated by 
the breaking index equations. 


MI. WEIGHTED LAPLACE EQUATION 


The Laplace equation is formulated by employing the 
continuity equation, substituting the properties of the 
velocity potential into the continuity equation. The 
formulation of the weighted continuity equation is carried 
out using a well-known method, namely by applying the 
principle of mass conservation to a control volume where 
fluid inflow-outflow occurs, 


3.1. Weighted Continuity Equation 


By employing the weighted Taylor series, a weighted 
continuity equation can be obtained. 


Vx = +y, = =O 0 (4) 

u is the horizontal water particle velocity, and w is the 
vertical water particle velocity. This equation represents the 
weighted continuity equation, where y, and y, are 
weighting coefficients. By using the weighted Taylor series, 
there is no longer truncation error or at least the truncation 


error has been greatly minimized 


Mathematically, equation (4) can be written as, 
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Both of these writings are somewhat inaccurate and will 
result in different solutions because each coefficient has a 


sis : : : : 4 
specific function. It can be said that yy is associated with = 


2 P +, dw 
and y, is associated with oe 
Z 


3.2. Weighted Laplace Equation 


In fluid flow, there is a scalar quantity called velocity 
potential, denoted as, 


Substituting the velocity potential property into (4), we 
obtain the weighted Laplace equation, 


a?e a?e 
Yx T2 +yz aa? =O 8 eee (5) 
This Laplace equation differs from the commonly used 
Laplace equation in water wave modeling, where there are 
actually weighting coefficients with yy = yz = 1. 


IV. SOLUTION TO THE WEIGHTED 
LAPLACE EQUATION 


4.1. Solution using the Separation of Variables Method 
The solution to (5) is carried out using the separation of 
variables method. In the variable separation method, it is 


assumed that the velocity potential is the product of three 
functions (Dean (1991)), namely 


o(x,z,t) = X(x)Z(z)Tt) ....(6) 
Here, X(x) : is a function of x only, Z(z) is a function of z 


only, and T(t) is a function of t only. Substituting into (5) 


and dividing the equation by (6), the following is obtained: 
Ye PX yz WZ 


— + — =l 
X(x)dx? Z(z)dz? 
This equation is fulfilled if 


Me TE na 
X (x) dx? 
Defined as k, = ee hence 
1 ax _ 2 
XG) ax? =k% aal) 


Using the same method, 


1 @Z_ 49 
zæ az2 k eens (8) 
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k 
Where k, = —. 
4 VYz 

. . 2 
k is the wave number where the wavelength is L = =< 
Hence, there are two wavelengths: the horizontal 


2 l 
wavelength, L, = = , and the vertical wavelength, L, = 
x 


= with different lengths. 


Equation (7), is offered for a solution 
X(x) = Acosk,x + B sin kx 
Equation (8) has a solution of 
Z(z) = Cekz? + De~ kz? 
Then, an assumption is made that the velocity potential is 
periodic with respect to time; hence, T(t) = sinot 

o= = s the angular frequency, where T is the wave 
period. Substituting X(x), Z(z) and T(t) ke (6), 

$(x,z,t) = (Acosk,x + B sink,,x) 

(Cez? + De-*#7)sinot — ....(9) 


The constants A, B, C, and D in equation (9) still need their 
specific forms to be determined. 


4.2. Working on the Kinematic Bottom Boundary Condition 


To obtain equations for the constants in the solution, the 
Kinematic Bottom Boundary Condition is applied at 
characteristic points where cosk,x = sink,x. At these 
characteristic points, the velocity potential equation 
becomes, 
o(x, z,t) = (A + B) cosk,.x 
(Cekz2 + De-*z?) sin ot 
Kinematic bottom boundary condition is, 
dh 


eS ae 


W_,» is the vertical water particle velocity at the sea bed at 
z=—h. 
U_p is the horizontal water particle velocity at the sea bed 
atz = —h. 

; : : dh . 
h is the water depth relative to the still water level, q IS the 


bottom slope, which has a negative value for waves moving 
from deeper to shallower waters. Utilizing the properties of 
velocity potential, 
do 
w(x,z,t) = Era —(A + B)k, cos kx 
Z 


(Ce¥z7 — De-*z7) sin ot 
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W_p = —(A + B)k, cos kx 
(Ce—*z" — Dekz") sin ot 


d 
u(x,z,t) = n = (A + B)k, sin kx 


(Cekz2 + De™*z7) sin ot 


u-n = (A + B)k, sink,x (Ce~*z" + De*z") sin ot 
The substitution of the Equations W_p and u_, to 
Kinematic Bottom Boundary Condition was performed at 
characteristic points as well as at (A + B) and and when 


sin(ot) is not equal to zero, yields the equation 


dh 
k,(Ce-*z" — Dekz) = k, (Ce—*z" + Dekz") ay 


dh dh 
(k -k —) Cenkeh = (k +k —) Dekeh 


* dx * dx 
Considering k, = £ dan k, = = 
g x V¥x Zz Vz 


1 1 dh ih 1 1 dh 
——— | Ce~*eh = +—— ] Dekh 
VV2 Jv 4x VYz JY- dx 


1 a 1 dh 
“Ty dx 
Pee i Vx x 
1 _ 1 dh 
Wz Jn 
Defined as, 
1, 1dh 
Vi dn” 
“=I Idh 
Wz Jn 
Hence 


C = De?kzhg 
Substituting to (7) 
p(x, z,t) = (A + B) cosk,x 


(De?kz wekz2 + De-*z7) sin ot 


(x, z, t) = (A + B)De*" cos kx 


(aekz@th) + e-kz@t) sin ot 


Defined as, 

aekz(Zth) +4 e Kz(zth) 
2 

aekz(zth) En e~Kz@EZ+h) 


2 


pE) = 


hC) = 
Where on a = 1, 


B(z) =coshk,(h+ z); B,(z) = sinhk,(h + z) 
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(x,z,t) = 2(A + B)De*2"B(z) cos kx sin ot 
Defined as A = 2A dan B = 2B, hence 
(x,z,t) = (A+ B)De*2"B(z) cosk,x sin ot 
Hutahaean (2022) shows that A = B 
(x,z,t) = 2 ADe*2"B(z) cosk,x sin ot 
Defined as G = ADe*z" 
o(x,z,t) = 2GB(z) cosk,.x sin ot 
The full equation, 

o(x, z,t) = GB(z) cosk,x sin ot + GB(z) 

sin k,x sin ot 


(10) 


At the characteristic point where cosk,x = sink, x, 
velocity potential equation becomes, 


b(x,z,t) = 2GB(z)cosk,xsinot —....... (11) 


vV. EQUATION FOR G 


The equation for G is formulated by integrating the 
Kinematic Free Surface Boundary Condition with respect to 
time t and is performed at characteristic points. The form of 
the Kinematic Free Surface Boundary Condition uses a 
weighted Taylor Series for two variables f = f (x, t), in this 
case, where the function f is the water elevation 7 = 
n(x, t). 


dn dn 
n(x + ôx, t + ôt) = n(x,t) + Ye2ôt r + YxÔX mr 
n(x + 6x,t + dt) — n(x, t) dn ôx dn 
ôt = Yez ae + Ye Se ay 
With ôt and ôx that are small, there obtained, 


Dn dn dn 


ae — + aa 
dt 2 ge * "nay 


Dn. : ; ; 
z is the total velocity of vertical surface water particle 


movement is denoted as w,, while u, represents the 
horizontal surface water particle velocity. Therefore, the 
equation for the Kinematic Free Surface Boundary 
Condition is. 

an 
"ax 
This equation is expressed as the water surface equation as 
follows, 


dn 
Wy = Yt,2 dt + YU 


n dn 
Yt,2 dt = Wy T Yxun dx 
Wy and u, are substituted by velocity potential (11), using 
the velocity potentials properties (u = — e w=— ze ; 
and by considering that k, = T dan k; = r 


Page | 37 


Hutahaean 


dn 1 dn 
Yt,2 ae —2Gk (Gam + V¥xB(n) a 


cosk,x sin(ot)  —....... (12) 


For the periodical function, 


2Gk (= A = tant 
E Bi) + VYxp n) Ax constan 
Integration (12) can be resolved by integrating sin(at). 


2G 1 d 
next) == (Thin + BO ae) 


k 
20 WYz 


cosk,.x cos(ot) 


For the periodical function, 


2Gk 


A= 2% (Lain) + Srb) 


Where A is wave amplitude. Hence, water surface elevation 


(13) 


equation is presented as, 


n(x,t) = A cosk,x cos(ot) (14 
dn : 
Pi —k,A sin k,x cos ot 


At the characteristic point, cosk,x =sink,x dan 
cos(ot) = sin(ot) 

dn kA 

dx 2 


Substituted to (13), obtaining a wave amplitude function 
equation as follows. 


a= alah) Wee) 


Considering k, = = and taking out £ (5) of the bracket, 


_ 2046 8) (1 60(8) ua 


Yt,27 VYz B (5) 2 


Based on the conservation law of wave number discussed in 


section (6), hence 6 (5) and 6, (5) are constant, where 


ps (3) 


= constant 


This constant value applies to all bodies of water, including 
both deep water and shallow water. In deep water where the 
bottom slope no longer has an effect, hence 


ps (3) 


(5) 


This condition is achieved when, 


A 
= tank (h+5) = 1 
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k(n+4) =6 
(1+7) =o 


Where tanh 8r ~ 1, @ is referred to as the deep water 
coefficient. The equation for the wave amplitude function 
becomes: 


2Gkcosh@n (tanhOn kA 
A= (Cae) (5 
Yt,20 VYz 2 (13) 
This equation can be expressed as Equation for G, 
_ OYt,2A 
G= TEN ZEZA hee (16) 
vyz 2 


VI. EQUATIONS OF CONSERVATION 


In its movement towards shallower waters, waves will 
undergo changes in its constants, namely G, k and A. The 
governing equations controlling these changes are the 
conservation equation of wave number and the conservation 
equation of energy. 


a. Conservation Equation of Wave Number 


In the method of separating variables, it is stated that Z(z) 
is a function of z only. For the velocity potential equation, 
Z(z)is expressed as: 


Z(z) = BP) 
As only function z, hence, 
dZ(z dk,(h+z 
dZ(z) = B,(z) Akak Ez) =0 
dx dx 
Hence 
dk,(h+z) _ i 
dx ~ 
Substituting k; = = and according to the formulation of G 
equation, therefore z = f is used. 
A 
ak(h+Ż 
2 e (17) 
dx 


This equation is the conservation equation of wave number. 
This equation shows that,, 


k (h + 2) = constant 
2 
The constant value is found in deep water, i.e., 
A 
tanhk(h+4) = 1 


Where tanh 0r ~ 1 and @ is the deep water coefficient, 
therefore, 


k(h+4) =6n (18) 


This equation is valid for all water depths, including deep 
water and shallow water. 


b. Conservation Equation of Energy 
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Substituting (10) into (5) using the properties of velocity 
potential, and working on the obtained equation with point 
characteristic properties, we get 

aG 
“ae 0 ...(19) 
Substituting (11) into (5) using the properties of velocity 
potential and considering (19) at characteristic points, we 
obtain 


GE +2 =0 ....(20) 
dx dx 


This is an energy conservation equation. 


Vil. DEEP WATER WAVE LENGTH 


Hutahaean (2022) formulates the relationship between 
wavelength and wave height in deep water. By employing 
weighting coefficients in this study, the relationship 
between wavelength and wave height in deep water is, 


Lg a 5 io (21) 


Ye3tanhon 


Lox = LU Vx 
Loz = L Syz 


Example calculation results of the wavelength for a wave 
amplitude A = 1.2 m, with various values of 


weighting coefficients, are presented in Table (2). 


Table.2 Wave length at wave amplitude of 1.2 m. 


E Ly Lz H 
(m) (m) Ly 
0.03 9.11 9.695 0.263 
0.032 9.161 9.839 0.262 
0.034 9.217 9.997 0.26 
0.036 9.279 10.172 0.259 
0.038 9.347 10.363 0.257 
0.04 9.421 10.573 0.255 
0.042 9.503 10.803 0.253 
0.044 9.592 11.055 0.25 


0.046 9.689 11.33 0.248 


In Table (2), the s values represent the weighting 
coefficients, with the corresponding values available in 
Table (1). 


Table (2) includes two wavelengths: the horizontal 
wavelength L, and the vertical wavelength L,. These 
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wavelengths exhibit a slight difference, although both are 
present. It is noteworthy that a larger € corresponds to a 
longer wavelength, resulting in a smaller wave steepness. 


Toffoli, A., Babanin, A., Onaroto, M., and Wased, T. 
(2010), determined that the critical wave steepness is 0.170, 
with a recommended upper limit of 0.200. To align with 
Toffoli et al.'s criteria and achieve a wave steepness closer 
to their recommendations, it is advisable to utilize a larger € 
value. 


VIII. WAVE TRANSFORMATION MODEL 


This section demonstrates that the potential velocity 
encompasses the shoaling-breaking phenomenon. Notably, 
this phenomenon persists irrespective of whether the 
original Laplace equation (as proposed by Hutahaean 
(2023b)) or the weighting coefficients yy = y, = 1 are 
employed. The focus of this section is solely on illustrating 
that the shoaling-breaking phenomenon remains unaffected 
by the weighting coefficients in the Laplace equation. 


Within this context, the weighting coefficient plays a crucial 
role in influencing the breaking characteristics. 
Specifically, it impacts key parameters such as breaker 
height H,, breaker depth h, and breaker length Ly. 


8.1. Shoaling Breaking Modeling 


The shoaling and breaking models are established based on 
conservation equations (17) and (20), incorporating 
equations governing the wave amplitude function (15) and 
the function G (16). The detailed formulation process is 
elucidated in Hutahaean (2023b). The equations governing 
the changes in wave constants k, A, and G as a wave travels 
from point x to x + 6x are essential components of this 
approach. 


dk _ 4k dh 


dx (4h+3A) dx 


dA_ G “(ao 


= “a cosh (67) (23) 


ax OYt,2 Ax 
Furthermore, 


dk 
kx+ôx a ky, + Ôx 
dA 
Ax+ôx = A, ae Ox a 


T 
ao = en Gy—5 (In ky4§x—In ky) 


To demonstrate the outcomes of the shoaling-breaking 
model, an analysis was conducted in a coastal zone 


characterized by a bottom slope of z = —0.02. The study 


utilized waves with a period of 8.0 seconds and a deep water 
wave amplitude of Ag = 1.2 meters. The deep water 
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coefficient 0 = 1.60, and a weighting coefficient of ¢ = 
0.042 was applied, with the specific values found in Table 


(1). The deep water wave number, kg, was calculated using 


equation (21), and the deep water wave height, họ = a = 
0 


A i 
oe Additionally, the deep water wave constant Gy was 


calculated using equation (16). The results of the shoaling- 
breaking model are visually presented in Figure (1), 
offering a comprehensive illustration of the dynamic 
behavior in the coastal zone under the specified conditions. 


12 
10 


H &L(m) 
ON FH WO 


0 2 4 6 8 
water depth h (m) 


wave height H wavelength L 


Fig.1: The outcome of shoaling-breaking modeling 


In Fig(1), it can be observed that the model is able to 
simulate shoaling and breaking phenomena effectively. 
Details regarding breaking characteristics, including 
breaker height Hp, breaker depth hp and breaker lengthL,, 
are discussed in another section. 


8.2 Refraction-Diffraction Analysis 


Utilizing the shoaling-breaking equations, the formulation 
of refraction-diffraction equations is established, as detailed 
in Hutahaean (2023b). The application of the refraction- 
diffraction model is demonstrated on bathymetry featuring 
a headland configuration (Fig. (2)), considering waves with 
a period of 8 seconds and a deep water wave amplitude of 
1.20 meters. The resulting 2-D contour image depicting 
wave height from the refraction-diffraction model is 
presented in Fig. (3), while Fig. (4) illustrates the 3-D wave 
height contour. 


/ $ 
s f as 


gt 


Lebel dL 


Fig.2: Batimetri tanjong contour 
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a deep water 


coefficient of 0 = 1.60 and an optimization coefficient of 
€ = 0.042, providing insights into wave behavior in the 
presence of bathymetric features and headlands. 


Fig.3: Wave height 2-D Contour 


Fig.4: wave height 3-D Contour 


In the analysis of refraction-diffraction results, a notable 
concentration of wave energy is observed at the center of 
the headland, coinciding with the point where breaking 
occurs. Both the shoaling-breaking model and the 
refraction-diffraction model incorporate equations 
involving the wave amplitude function, wave number 
conservation, and energy conservation to collectively 
simulate the processes of shoaling and breaking. These 
equations include weighting coefficients. 


To determine suitable values for these weighting 
coefficients, an examination will be conducted using 
breaker index equations, as elaborated upon in the 
subsequent section.. 
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IX. BREAKING INDEXES EQUATIONS. 


In the preceding section, it has been demonstrated that the 
developed equations can effectively simulate shoaling and 
breaking. In this section, breaking indexes equations are 
formulated, and the role of corresponding weighting 
coefficients is studied. 


9.1.Equation breaking. 


The breaking equation is present in both the wave amplitude 
function Equation (15) and the wave constant Equation G 
(16). Breaking occurs when, 
tanh6x kA _ 
We 2 
Based on Equation (24), breaker indexes are formulated, 
including the breaker length index, breaker depth index, and 
breaker height index. 


Oo Ss (24) 


9.2. Breaker Length Index Equation 


Equation (24) can be expressed as, 


kA _ tanh ém 
2 VYz 
Considering k, = = or k =k,../y, and ky, = = hence 


Hp _ 2tanh 6n/yx 
lig ae OOO CS 


H, represents the breaking wave height, while Lyp is the 
horizontal wavelength at the breaker point. In Equation 
(25), there are weighting coefficients, namely yy and yz, 
indicating the influence of weighting coefficients on 
breaking characteristics in this equation. In Table (3), the 


: Hp . : 
calculation of aa is presented for various values of € and 
b,x 


6, where 6, = 1.60, 0, = 1.70, 03 = 1.80 and 6, = 1.90. 
The values of tanh 8r include, 


tanh 1.67 = 0.999914 
tanh 1.77 = 0.999954 
tanh 1.87 = 0.999975 
tanh 1.97 = 0.999987 
All of them ~ 1. 


Table.3 The value of — on different € and 0 
b,x 


E A 0, 03 0, 

0.03 | 0565] 0.565) 0.565| 0.565 
0.032 | 0.555 | 0555] 0.555] 0.555 
0.034 | 0.545 | 0.545] 0.545] 0.545 
0.036 | 0.535 | 0535] 0.535] 0.535 
0.038 | 0.524 | 0.524| 0.524] 0.524 
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0.04 0.513 0.513 0.513 0.513 
0.042 0.502 | 0.502 0.502 0.502 
0.044 0.49 0.49 0.49 0.49 
0.046 0.479 | 0.479 0.479 0.479 


In Table (3), it can be observed that as the value of € 

: H : 

increases, the value of z decreases. Meanwhile, 
b,x 


‘ at H 
concerning the deep water coefficient 0, the value of n 
b,x 


remains constant, as this is only considered up to the third 
decimal place. If examined up to the fifth decimal place, the 


H : : : ; 
aE value increases with an increase in the 0 value. 
b,x 


9.3. Breaker Depth Index Equation. 
Defined as 


t nk(n+5) = k(n+5) 
an 7 =a 7 


Substituting to (24), 


ak(h+5) ka 


0 
VYz 2 
A VVzkA 
k(n Z) - =0 
a t3 7 
akA kA 


A 
a -— (yy - a); =? 
k(n+4) = Or hence , 


Substituting to Equation a, 


E tanh 67 
— on 


At breaker point and on the sinusoidal wave, A = z applies 


tanh 0r tanh 2) A 


On -(Ve- Ox /2h 


H 
tanh 67 — (Ony; — tanh 0r) =0 
b 


Hb _ 4 tanh 0r 
hp  (Onvyyz—-tanh 0r) 


This equation represents the breaker depth index, where h, 
is the breaker depth. In equation (26), there is a weighting 
coefficient denoted as y,. Table (4) presents the results of 
the calculation of equation (26) for various values of @ and 
several values of £. In the table, 04 = 1.60, 862 = 1.70, 63 = 
1.80 and 6, = 1.90. 
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Table.4 the value of 2 on several values of € and 0 
b 


E A 0, 0; 0, 

0.03 | 0.869] 0.808] 0.754] 0.708 
0.032 | 0.853 | 0.793] 0.741 | 0.695 
0.034 | 0.836 | 0.777] 0.726| 0.682 
0.036 | 0.819 | 0.762] 0.712] 0.668 
0.038 | 0.801 | 0.745] 0.697] 0.654 

0.04 | 0.783 | 0.728 | 0.681 | 0.639 
0.042 | 0.764| 0.711] 0.665] 0.625 
0.044 | 0.745 | 0.694] 0.649 0.61 
0.046 | 0.726 | 0.676 | 0.633 | 0.594 


The results in Table (4) indicate that the value of %2 
b 


decreases with an increase in € and also decreases with an 

p ; . H 

increase in 8. For a constant value of H}, the decrease in = 
b 


with an increase in @ suggests a deeper breaker depth. 


According to Mc. Cowan's criteria (1894), where = = 
b 


0.78, and considering the calculation of the breaker height 
h, in subsection (9.5), resulting in € = 0.038 — 0.046, the 
appropriate 0 value is therefore 1.60. 


9.4 Breaker depth-length index. 
Equation (25) is written as an equation for H, and 


substituted into (26), resulting in the breaker depth-length 
index equation, denoted as Equation (27). 


I» E — et) NA 27) 


Lb,x 


nelh oe ; 
Condition z > 1 represents deep water conditions. Breaking 


due to bathymetry occurs in shallow water where t< 1 
hence, (27) must be less than 1. Table (5) illustrates the 
values of “ for = 1.60 and 8 = 2.30 

b,x 


: . h 
Table.5 Provides an overview of the values of 
b 


concerning the deep water coefficient 0. 


E 6 = 1.60 6 = 2.30 
hp H, hp H, 
Lox hy Lyx h, 


0.03 0.65 0.869 0.996 0.567 
0.032 0.651 0.853 0.997 0.557 
0.034 0.652 0.836 0.997 0.547 
0.036 0.653 0.819 0.998 0.536 
0.038 0.655 0.801 0.998 0.525 
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0.04 0.656 0.783 0.999 0.514 
0.042 0.657 0.764 0.999 0.502 
0.044 0.658 0.745 1 0.491 
0.046 0.659 0.726 1 0.479 


From Table (5), it is evident that as the value of € increases, 
h f i 

the value of an , also increases, and as @ increases, the value 
b,x 


h : h 
of —2 also increases. At 0 = 2.30, the value of —2 
Lyx Lyx 


approaches 1. It can be concluded that the value of 0 must 
be less than 2.30. 


9.5. Breaker Height Index Equation, 2 
0 


The breaker height index is the ratio of the breaker height to 
the wave height in deep water, denoted as > The wave 
energy for a single wavelength is given by the equation: 

E = cgpgH?L 

where cp is the energy coefficient. In the linear wave theory 
(Dean (1991)), cg = =. In this context, the value of cz is not 


influential as it cancels out in the energy conservation 
equation. By equating the wave energy at the breaker point 
to the wave energy in deep water, we obtain Equation (28): 


HÈL, = Hlo - |. -xpeeee (28) 
Equation (25) can be expressed as an equation for Lp and 


substituted into (28). Meanwhile, Lọ is substituted with 
(21), resulting in Equation (29): 


1 
Hy _ (2 Ya (ve2+ 23) I (29) 
Ho Yt,3 VYz E 


In this equation, there are no parameters for the deep water 
coefficient 0 and wave period. Therefore, the breaker height 
is solely determined by the deep water wave height Họ and 
the weighting coefficient. 


To illustrate the influence of the weighting coefficient on 

the breaker height index “ and the breaker height Hp, the 
0 

results are presented in Table (6). The calculated breaker 


height is compared with the breaker height from Komar, 
Paul D., and Gaughan, Michael K. (1968): 


H, = 0.39 g /S(THo) /5 m. 


In this equation, the wave period T, is set to 8 seconds, and 
two different deep water wave heights Hy of 2.00 m and 
2.40 m are used. In Table (6), Hj, —29 is the result of Equation 
(29) multiplied by Ho. 


In Table (6), an observable trend is noted: as € increases, the 
breaker height decreases. Following the Komar-Gaughan 
equation, a breaker height Hy = 2.00 m is achieved at ¢ = 
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0.038, while for a higher initial wave height Hy = 2.40 m, 
the corresponding ¢ = 0.046. However, it is challenging to 
precisely determine e from these findings. For a wave height 
of 2.40 m and a wave period of 8 seconds, identified as the 
maximum in that period, a pragmatic approach suggests 
using € = 0.042 for wave heights below the maximum. 


Table.6 Presents the calculated breaker height Hp. 
E H, Ho = 2.0m Ho = 2.4m 
Ho Hp-29 | Hp- | Hp-29 | Hp-x 
0.03 1.267 | 2.534 | 2.463 | 3.041 | 2.85 
0.032 | 1.259 | 2.518 | 2.463 | 3.021 | 2.85 
0.034 | 1.25 2.501 | 2.463 | 3.001 | 2.85 
0.036 | 1.241 | 2.483 | 2.463 | 2.979 | 2.85 
0.038 | 1.232 | 2.464 | 2.463 | 2.957 | 2.85 
0.04 1.222 | 2.444 | 2.463 | 2.933 | 2.85 
0.042 | 1.212 | 2.424 | 2.463 | 2.909 | 2.85 
0.044 | 1.202 | 2.403 | 2.463 | 2.884 | 2.85 
0.046 | 1.191 | 2.382 | 2.463 | 2.858 | 2.85 


The study of breaker indexes underscores the substantial 
impact of weighting coefficients on breaking parameters. 
Consequently, employing the weighted Laplace equation in 
modeling wave dynamics is anticipated to yield more 
accurate and reliable results for breaking parameters. 


The results of this study also indicate that the appropriate 
deep water coefficient is 0 = 1.60, with a value of € = 
0.042, along with weighting coefficients: 


E Yt2 Yt,3 Yx Yz 
0.042 | 1.99626 | 3.1097 | 0.97831 1.26448 


X. CONCLUSION 


The Weighted Taylor series refers to a truncated Taylor 
series that is transformed into a first-order series with 
significantly reduced truncation errors, and in some cases, 
without any errors. Utilizing the weighted Taylor series for 
formulating the continuity equation and Laplace equation 
results in the creation of the weighted continuity equation 
and weighted Laplace equation, characterized by minimal 
truncation errors. Consequently, constructing a water wave 
model using the weighted Laplace equation enhances the 
model's accuracy. 


The pivotal role of weighting coefficients cannot be 
overstated, as they play a crucial role in determining wave 
characteristics such as wavelength and breaking parameters, 
including breaker height, depth, and length. This 


www.ijaers.com 


International Journal of Advanced Engineering Research and Science, 11(1)-2024 


underscores the importance of utilizing the weighted 
Laplace equation in water wave modeling. Furthermore, 
since water wave modeling is inherently based on velocity 
potential, employing the weighted Laplace equation is 
deemed more effective. 


In situations where the latest data from physical model 
research is available and reliable, the values of weighting 
coefficients can be easily adjusted to fine-tune the model. 
This adaptability enhances the applicability and accuracy of 
the water wave model, making it a valuable tool in studying 
water wave behavior. 
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